Responses of quark condensates to the chemical potential 
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The responses of quark condensates to the chemical potential, as a function of temperature T and 
chemical potential fi, are calculated within the Nambu-Jona-Lasinio (NJL) model. We compare our 
results with those from the recent lattice QCD simulations [QCD-TARO Collaboration, Nucl. Phys. 
B (Proc. Suppl.) 106, 462 (2002)]. The NJL model and lattice calculations show qualitatively 
similar behavior, and they will be complimentary ways to study hadrons at finite density. The 
behavior above T c requires more elaborated analyses. 
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PACS numbers: 12.39.-x 

The variation of quark condensates in medium plays 
a key role to understand the behavior of hadron masses 
and chiral symmetry restoration ^ . Recently, we cal- 
culated for the first time the second order response of the 
quark condensate to the chemical potential d ''Jup at fi = 
using lattice QCD || following the method in M. It 
was found that the response is negative both below and 
above T c . One of interesting results is that the response 
to the isoscalar chemical potential (/j,g = fi u — fid) is al- 
most the same as that to the isovector chemical potential 
(pv = Mm = ~ fJ-d) , where fi u (fid) is the u (d) quark chem- 
ical potential. It would be interesting if we can check this 
result within effective models of QCD. In this work we 
present an Nambu-Jona-Lasinio (NJL) model H |(| cal- 
culation of **£ 9 ' and d at fi = and compare our 
results with those from the lattice QCD simulations ||. 

An SU(2) NJL model Lagrangian will be enough for 
that purpose. However, we are also interested in the 
responses of the s-quark condensate and a comparison 
with those of the u,d-quark condensates will be useful for 
future studies in the lattice calculations. Thus we used an 
SU(3) NJL model in our calculations and found that the 
effects of the flavor mixing are negligible. We present only 
the results for the u,d-quark condensates in this paper, 
and we will make the comparison in a forthcoming paper. 

First, let us consider an SU(3) NJL model Lagrangian 
0: 

1 8 

L = q(ij ■ d-m)q + -g s ^2[(qX a q) 2 + {qi\ a j 5 q) 2 ] 



+ g D [det qi(l - -f 5 )q 3 + h.c 



m d = 5.5 MeV, 



135.7 MeV, 



(2) 



where A is the momentum cut-off. The third term in 
Eq.(|l]) is a reflection of the axial anomaly, and causes 
a mixing in flavors. For example, the constituent quark 
masses are given as follows. 



M u = m u -2g s (uu) - 2g D (dd)(ss), 
M d = m d - 2g s (dd) - 2g D (uu)(ss), 
M s = m s - 2g s {ss) - 2g D (uu)(dd), 



(3) 



where (•) means the statistical average. In this work we 
concentrate mostly on the Case II in 0, where only go 
has a temperature dependence 



g D (T) =g D (T = 0) exp[- 



-(T/To) 2 ] 



(4) 



while other coupling constants and the cut-off are inde- 
pendent of T and chemical potential (or density). Here, 
we set To = 0.1 GeV taking into account the restoration 
of Ua(X) symmetry as in [Q. 

In the mean-field approximation the above Lagrangian 
leads to the following gap equation [pj: 



Q t =2N c J2 



-Mj 



w. 



(QiQi), 



(5) 



where the index i denotes the u, d, and s quarks. N c 
is the number of colors and Mi is the constituent quark 



mass, and E ip 

nip and n 



where 



VW+pt. f(E ip ) 



1 



ip 



are the distribution functions of the 



(1) ith quark and antiquark, respectively, i.e. 



where A a are the Gell-Mann matrices and m is a mass 
matrix for current quarks, m=diag(m u , md, rn s ). We 
take the following parameters in : 

A = 631.4 MeV, 5s A 2 = 3.67, g D A 5 = -9.29 



'Deceased. 



n 



1 



1 + cxp {{E ip -fii)/Ty 
1 

'" " 1 + cxp {{E w + fi t )/T) ' 



(6) 



The l.h.s. of Eq.(g) is a function of (uu), (dd), (ss), fii, 
and T. Then, we obtain the first order response of quark 



condensates 



d{uu) 
dp 



d(dd) 
dp 



and 



d(ss) 
dp 



by differentiating both 



2 



sides with respect to /i at a fixed T, i.e., we solve the 
following equations 



dQ u ((uu), (dd), ■ 




d(uu) 


dfi 




dfi 


dQ d ((uu), (dd), 




d(dd) 


d/i 




dfi 


dQ s ((uu), (dd), 


(ss),n s ) 


d(ss) 


dfi 




dfi 



(7) 



Here, we consider two types of u following the notation 
in the lattice QCD simulations ||, ||| . One is the isoscalar 
chemical potential fis = fiu + Mrf, and the other is the 
isovector one fly = fi u - fid- 
Let us comment on the definition of chemical potential 
used in the lattice calculations and ours. There will be 
no difference between the two definitions if we choose the 
same u- and d-quark mass. In the NJL model calculations 
we chose rather general definitions for the isoscalar and 
isovector chemical potentials, i.e., fis = Mu + an d 
= fi u — fid to show the quark mass dependence by taking 
different u- and d-quark masses. 

Using Eq. (0) one can show that both d j^} and d j^ 
are zero at £t=0. This is the same as the lattice QCD 
calculation ||. At finite chemical potential, the absolute 
value of the quark condensate decreases with increasing 
chemical potential and the variation is proportional to 
the chemical potential in the present NJL model J§| . 

Next, consider the second order response d d ^ ■ This is 
obtained by differentiating each equation in Eq. ([?]) with 
respect to fis (or fiy) again. In Fig [l] we show d d ^£^ 



for the u and d quark condensates at fi = 0. Within 
the present NJL model the second order response to the 
isoscalar chemical potential fig is the same as that to 
the isovector chemical potential fly . This is consistent 
with the result in ||. In fact, there are a few different 
terms between the isoscalar and the isovector cases in the 
lattice QCD simulations. However, it is found that the 
contribution of those terms is negligible as shown in the 
figure 3 of Ref.[§. 

In the lattice calculations the same u- and d-quark 
mass was taken ||. In this paper, however, we would 
like to show the quark mass dependence of the response in 
the figure. Our main interest is to compare the responses 
below and above the Mott temperature depending on the 
quark mass. The present results may be useful for future 
lattice QCD simulations which use different u- and d- 
quark masses. 

We take m u — 4 MeV and rrid — 7 MeV as well as m u 
— m d = 5.5 MeV to show the quark mass dependence of 
d ^Sp. I , Although the cut-off and the coupling constants 
should be modified, we use the same coupling constants 



and the cut-off for both cases and study the behavior of 
d j2 ■ Below the pion Mott temperature T m ^ (« 148 
MeV in this work), the response of the u (m M = 4 MeV) 
quark condensate is larger than that of the d (md = 7 
MeV) quark condensate, while above the Mott tempera- 
ture the behavior of the response is opposite (Here, T m?r 
is determined as a temperature at which the sum of the u 
and d constituent quark masses equals to the pion mass, 
i.e., M u + Md = m^). This behavior results from the 
quark mass dependence of the first order response at fi- 
nite chemical potential shown in Fig. ^, where we take fi u 
= fid — 0.02 and 0.04 GeV, respectively. That behavior 
in Fig. H is not altered even in the Case I 0, where all 
the couplings and the cut-off are independent of temper- 
ature and chemical potential. For example, in the Case 
I and at fi — 0.02 GeV, the pion Mott temperature is 
about 203 MeV and the extremum is located at around 
186 MeV. 

In Fig. |l| we also present the relative variation 

d du? I (99) at ^— f° r comparison. Now, d }^ / (<?g) has 
the extremum at the Mott temperature and the (absolute 
value of) variation decreases with increasing the quark 
mass. We find similar behavior in the lattice data, al- 
though (qq) is very small over T c in this case. 

As we know, we can not detect a quark condensate 
itself in experiments. However, variations of quark con- 
densates, such as shown in Fig. ||, are necessary to pre- 
dict responses of hadron masses to the chemical potential, 
e.g., in ||. A more detailed analysis on the extremum 
point and its effects on hadron masses is in progress. 

In summary, we have calculated the first and second 
order responses of quark condensates to the chemical po- 
tential within the NJL model and found that they are 
consistent with those from the recent lattice QCD simu- 
lations. 
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FIG. 1: The second order responses of quark condensates at fi = 0: J-™' (left) and /(gg) (right). Here, means both 
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FIG. 2: The first order responses of the u and d quark con- 
densates at finite chemical potential, where m u — 4 MeV and 
m d = 7 MeV. 



